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ABSTRACT 

We introduce the Fast Holographic Deconvolution method for analyzing interferometric radio data. 
Our new method is an extension of A-projection/software-holography/forward modeling analysis tech- 
niques and shares their precision deconvolution and widefield polarimetry, while being significantly 
faster than current implementations that use full direction-dependent antenna gains. Using data from 
the MWA 32 antenna prototype, we demonstrate the effectiveness and precision of our new algorithm. 
Fast Holographic Deconvolution may be particularly important for upcoming 21 cm cosmology ob- 
servations of the Epoch of Reionization and Dark Energy where foreground subtraction is intimately 
related to the precision of the data reduction. 



1. INTRODUCTION 

Interferometric radio arrays are some of the most scien- 
tifically productive tools of astrophysics. In particular, 
by choosing the appropriate antenna configuration the 
angular sensitivity of an interferometer can be tuned to 
the scientific question at hand. However, this flexibility 
comes at a the cost of increased demands on the precision 
of the data analysis. 

Radio interferometry uses an array of antennas to sam- 
ple the electric field from cosmic sources. Each antenna 
collects the electric field across its surface, which is then 
cross-correlated against all of the other antenna measure- 
ments to generate pairwise visibilities. Each visibility is a 
measurement of an angular Fourier mode of the sky (for 
an introduction to interferometry see Rohlfs & Wilso3 
[200 g) . The sparse sampling of angular Fourier modes 
allows a unique combination resolution and angular sen- 
sitivity tailored to the science case, at the expense of a 
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complicated point spread function. Additionally, each 
antenna integrates the electric field across the antenna 
surface — fundamentally erasing information below the 
antenna scale and imprinting a direction-dependent gain 
and polar ization response t o the signal (see chapter 3 
of iMorales fc Wvithg 120101 for a conceptual introduc- 
tion to dir ection-depende nt antenna effects). Recent ad- 
vances by iBhatnagar et al. (2008|), IMorales &: MateieH 
(|2009t ). and lMvers et al.l (,2003) — variously referred to as 
"A-projection" or "software holography" — have enabled 
these direction-dependent gain and polarization effects 
to be included in the deconvolution process. While the 
resulting analyses are quite precise, the y are slow com - 
pared to traditional CLEAN methods (IHogbomI [l974D . 
particularly if all the antennas do not share the same 
gain pattern. 

In this article we introduce the Fast Holographic De- 
convolution technique which increases the speed of A- 
projection/software-holography deconvolution. In iJ5]we 
introduce the Fast Holographic Deconvolution algorithm 
and compare it to commonly used radio deconvolution 
algorithms. We then describe our particular implemen- 
tation of Fast Holographic Deconvolution in fJSl and ap- 
ply the algorithm to data from the 32 antenna MWA 
prototype in fjll We conclude with a discussion of future 
work in fJS] 

2. ALGORITHM 

There are many variations of radio deconvolu tion in the 
literature, including CLE AN (jH ogbom' HJTI. Cotton- 
Schw a b facete d CLE AN (jSchV ab 1984: iCotton et all 
I2OOI ICottonI [2005I) . multi-scale CLEAN, MEMs 
A-projection based C LEAN (Bhatnagar et al.l |2008|; 
ICarozzi k WoanI I2009D. peeling (iMi tchcU ct al. 200^ 



forwa rd-modelling ( Bernardi et al.l 1 2011; .Pindor et aLl 
l2011f ). measuremen t equatio n (Smirnovl l2011a[ ). and com- 
pressive sampling ( jLi et al.l [201 ll) . to name just a few. 
Very crudely, the most popular styles of deconvolution 
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can be divided into two families: basic CLEAN im- 
plementations that subtract model visibilities from the 
data, and peeling, A-projection, and forward modeling 
approaches that subtract improved visibility estimates 
from the data with direction-dependent antenna calibra- 
tion. While this overly broad classification does not do 
justice to the many important advances in deconvolu- 
tion, this simplification helps place the advances of Fast 
Holographic Deconvolution in context. 

2.1. CLEAN 

The basic CLEAN alg orithm loop — omitting for clarity 
widefield w-projection (jCornwell et al.l I2008D . faceting, 
major-minor cycles and other innovations — can be ap- 
proximated with the following equation to determine the 
next component of the model: 

7Maxi F'^(6',u)S(u,v) 



G(v, v)vrf - 




(1) 



where we use the linear algebra notation A(a, b; x) to 
signify a transformation from coordinate vector b to co- 
ordinate vector a that depends on parameters x, and 
vector u is the uv plane and v represents a list of visibil- 
ities. Describing Equation [T] in words, the data visibili- 
ties Vd are calibrated (G) and the model visibilities v„j 
for all the previously found components are subtracted. 
The terms in parentheses are then the residual visibili- 
ties which are gridded to the uv plane (S) and Fourier 
Transformed (F ) to form the residual dirty map. The 
peak of the map is selected and the loop gain 7 is applied 
to determine the next component of the model and the 
loop is repeated. 

The model visibilities are calculated by Fourier 
transforming the model source components I {0m) to the 
uv plane and taking the (5-function of the correlation 
function at the baseline separations, as shown in this 
annotated and expanded version 

gridding calibrated data 

7MaxJ F^((?,u)S(^| G(^^^^ - 



model visibilities 



^<5(v,u)F(u,0)/(a,„) 



(2) 



This is of course a gross simplification of the algorithms 
in common use, but it captures the three key features of 
most CLEAN implementations: the model visibilities are 
the component correlations at the baseline separation, 
source components are at pixel centers with a magnitude 
determined by a CLEAN gain 7 that is less than one, 
and the model is subtracted from the data as visibilities. 

2.2. A-projection 

A-projection makes tw o significant chan ges to the 
basic CLEAN algorithm (jBhatnagar et al.l 12008). A- 
projection uses the holographic beam pattern (B) to grid 
residual visibilities to the uv plane, effectively applying 



direction-dependent calibration, and it uses the beam 
pattern (B) to calculate more accurate model visibili- 
ties. 



jMnxi F^ (6l,u)B (u,v 



Vd-^B(v,u)F(u,0)/(^. 




(3) 



In words, the model components I{9m) are transformed 
to the uv plane which is integrated by the known an- 
tenna response B(v,u) to form an improved estimate of 
the visibilities. Many array simulators such as MAPS 
(MIT Array Performance SimulatoiQ) operate by per- 
forming the uv integration described by B(v, u). This 
improved model visibility is subtracted from the mea- 
sured data, and the residual visibilities are gridded using 
the A-projection/software holography kernels B to form 
a holographic dirty map. The peak of the image is used 
to produce the next model component, and the process 
is repeated. 

In some situations it is faster to apply direction depen- 
dent gains (G(v, v; a, 6m)) to each model visibility rather 
than integrate the uv plane with B(v, u), giving 



7MaxCentroid< F^ (6',u)B (u,v) 



l^vd - G^(v, v; a, 0m)S{^, u)F(u, 0)/(0„O j | ■ (4) 

Mathematically these are identical as B(v, u) = 
G(v, v; a, 0m)<J(v, u), but they can have very differ- 
ent speeds depending on context. Using the direc- 
tion dependent gain is sometimes referred to as peeling 
(jMitchell et al.l 12008) . This is also similar to the mea- 
surement equation ap proach to directio n-dependent cal- 
ibration developed bv iSmirnovl ()2011bl O). 

2.3. Forward Modelling 

For ward modeling (jBernardi et al.|[2011l : IPindor et al.l 
120111) is mathematically similar to A-projection and peel- 
ing, but subtracts holographic model and data images 
instead of the model and data visibilities: 



7 MaxCentroid'l 



,u)B (u,v)vd 



(u,v)B(v,u)^F(u,0)/(0„)i. (5) 



The data is gridded and imaged [shown in square brack- 
ets] to form a holographic image (first line), and the 
model is converted to visibilities then gridded and im- 
aged to form a holographic model image (second line), 
which are then subtracted. For traditional dirty im- 
ages image subtraction would cause a significant degra- 
dation in deconvolution fidelity. However, because the 
antenna beam pattern (B) is used to grid the data, 
holographic dirty maps contain all the information that 

^ |http: //www. haystack. mlt . edu/ast/arrays/maps/lndex .html] 
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was in the visibilities and subtracting images or visibil- 
ities works equally well (assu ming calibration B is cor- 
rect, |Morales^^ateiel?'200^. This lossless information 
property was proved by Tcgmark (1997) and is widely 
used for reducing CMB data. One advantage of sub- 
tracting images is that the data only needs to be gridded 
and imaged once [square brakets] . Forward modeling also 
usually centroids the position of the source components 
as opposed to using the pixel centers, though component 
centroiding can be added to any of the deconvolution 
approaches. 

2.4. Fast Holographic Deconvolution 

The major change behind Fast Holographic Deconvo- 
lution is to introduce the Holographic Mapping function 

H(u, u) = B (u,v)B(v, u). In all of the holographic 
deconvolution algorithms (A-projection Eq. [3l peeling 
with holographic gridding Eq. |4l and Forward Model- 
ing Eq. [5]), the majority of the time is spent integrating 
regions of the model uv plane to form the model visibili- 
ties (B(v, u)), and then gridding the visibilities with the 

~ T 

holographic antenna beams (B (u, v)). Instead of re- 
peating these operations for every model component and 
visibility, we combine them into a Holographic Mapping 
function. Conceptually the Holographic Mapping func- 
tion H(u, u) records how a uv location in the model is 
mapped to the holographic uv map, independent of the 
sky model used and including all of the baseline sam- 
pling and direction-dependent beam effects. Effectively 
we pre-compute the process of forming and gridding visi- 
bilities and use this repeatedly as we refine the sky model. 
The Holographic Mapping function allows us to quickly 
convert a uv model into an accurate estimate of the holo- 
graphic uv map as seen by the instrument, including all 
of the direction-dependent array beam, antenna beam, 
and polarization effects. The resulting Fast Holographic 
Deconvolution algorithm is: 



7 MaxCentroid-l 



~ T 

,u)B (u,v)vd 



?,u)H(u,u)^F(u,0)/(0„ 



(6) 



In words, the source components are used to build a 
uv model (X^m '^("i ^)^(^m)), the Holographic Mapping 
function converts this to a model holographic uv map 
including all of the instrumental and gridding effects, 
a Fourier transform creates a model holographic dirty 
image (lower line) which is subtracted from the dirty 
holographic map of the data [square brackets]. Mathe- 
matically FHD is equivalent to the A-projection, peeling 
and Forward Modeling deconvolution algorithms^ but 
is more efficient since the effect of individual direction- 
dependent antenna gains is pre-computed outside of the 
deconvolution loop. 

The speed of Fast Holographic Deconvolution is due 
to the sparseness of the Holographic Mapping function. 
Formally H(u, u) is a very large matrix as it can map any 
model uv point to any point in the holographic uv map. 

^ Implicit gridding of the model uv plane can introduce small 
errors far from the field center. 



However, every visibility is the integral of a very limited 
region of the uv plane — about twice the diameter of an 
antenna. The gridding function is similarly compact. So 
in practice a location in the model uv plane is mapped 
only to nearby locations in the holographic uv plane. 
This combines with the proximity of nearby frequency 
and time samples to make the Holographic Mapping ma- 
trix H(u, u) very sparse and amenable to sparse matrix 
computational methods. Determining the sparse Holo- 
graphic Mapping function takes about ten times longer 
than one integration/gridding cycle (Table [Hand associ- 
ated text). But after H(u, u) has been calculated, we can 
take any model uv contribution and very quickly map it 
into a holographic dirty map with high precision. 

The power of the FHD technique comes from combin- 
ing image based subtraction and the Holographic Map- 
ping function. This enables us to produce a precision 
deconvolution method, complete with direction and po- 
larization effects, that is faster than current implemen- 
tations. 

3. IMPLEMENTATION 

In this section we detail our implementation of Fast 
Holographic Deconvolution. The FHD algorithm has 
three key steps: creating a holographic dirty image of 
the data, calculating the Holographic Mapping function 
H(u, u), and iterating through the deconvolution loop 
(Equation [5]). 

3.1. The Holographic Dirty Image 

We begin Fast Holographic Deconvolution by con- 
structing a dirty holographic image of the calibrated 
data for each of the 4 instrumental polarizations. We 
assume a dual linear reception element (labeled X and 
Y) , and create separate holographic dirty maps from the 
XX, YY, XY and YX visibilities. The maps are cre- 
ated by gridding each visibility with a kernel based on 
the individual holographic antenna maps for each polar- 
ization, and Fourier transforming the result. It is the 
holographic gridding kernel that enables the correction 
of the direction and polarization dependent response of 
each antenna. Our implementation supports arbitrary 
antenna kernels and any desired bandwidth synthesis. 
Since we use image based subtraction (Equation [S]), we 
only need to grid the visibilities to the uv plane and EFT 
once for each instrumental polarization. 

3.2. The Holographic Mapping Function 

The Holographic Mapping function H(u, u) must be 
pre-computed for every unique observation and for each 
instrumental polarization. Once computed, it is stored 
in the row-indexed sparse storage mode from Numerical 
Recipes in C and written to disk. Despite being very 
sparse, this is still typically a large file — a few gigabytes 
for our example in §4. However, the size of the matrix 
scales with uv coverage rather than the number of base- 
lines or the length of integration, so we believe this ap- 
proach will remain practical even for future large arrays. 
Finally, we note that the choice of pixelization of the uv 
plane is quite flexible. For example, the dimensions and 
resolution of the input and output uv planes do not have 
to match, or one could be a standard grid and the other 
non-rectilinear such as a uv analog of Healpix. 
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3.3. Deconvolution 

Once the Holographic Mapping Function has been 
built we begin the deconvolution loop described math- 
ematically by Equation [S] In words, the steps within 
each iteration are: 

1. We begin by applying the Holographic Mapping 
function to the model sky (in uv space) for each 
instrumental polarization, and take the FFT to cre- 
ate four model dirty holographic images. 

2. These are subtracted from the corresponding dirty 
holographic images of the data. 

3. The residual images are converted to Stokes I, Q, 
U, and V using the Jones matrices and an average 
beam model for each instrumental polarization. 

4. The Stokes images are median filtered with a small 
box size to highlight point sources and masked to 
exclude regions of sky with low instrumental re- 
sponse. (Median filtering is equivalent to down 
weighting short baselines, but is faster for holo- 
graphic algorithms.) 

5. The positive peak in the filtered Stokes I image is 
identified. A centroid is fit around the peak and 
the gain factor 7 is used to generate a new source 
component. (Q, U, and V amplitudes of the com- 
ponent may optionally also be fit.) 

6. The new component is used to update the uv plane 
model sky in all four instrumental polarizations. 

When deconvolution is finished, we condense all com- 
ponents within a threshold radius of half a pixel to form 
a list of discrete point source candidates. A robust source 
list may be constructed by applying standard statistical 
tests to the list of candidates, such as imposing a mini- 
mum signal to noise threshold. Extended sources which 
have more than a negligible number of close components 
outside of this threshold may be either left in component 
form or fit with a model profile. 

3.4. Comparison of computational efficiency 

Comparing the computational efficiency of very dif- 
ferent algorithms is always difficult due to the details 
of implementation and computer hardware. Instead of 
counting the formal number of operations, which tends 
to ignore important effects such as memory bandwidth, 
we have created matched implementations of both Fast 
Holographic Deconvolution and A-projection. The IDL 
programming language was used for both, and no par- 
allelization or optimizations were used (many optimiza- 
tions work equally well for both). Once we finish algo- 
rithm development we intend to translate our FHD im- 
plementation to C or CUDA and parallelize as necessary. 

The computational cost of deconvolving N components 
with A-projection can be represented by 



TA ~ N (rniodcl + Tgrid + Tfit) 



(7) 



where Tmodci is the time spent constructing model visibil- 

ities (B (u, v)), Tgiid is the time spent gridding the visi- 
bilities using the individual direction-dependent antenna 



gains to the uv plane (B(v, u)), and Tat is the time spent 
determining the next model source component. Forward 
modeling (Equation [5]) is computationally similar in cost 
to A-projection because it requires constructing model 
visibilities and gridding them back to the uv plane for 
each model component. 

In Fast Holographic Deconvolution there is an up-front 
cost for calculating the Holographic Mapping Function, 
followed by more efficient determination of each model 
component. The computational cost of FHD is repre- 
sented by 



Thd 



S (t„ 



Ddol 



+ Tgrid) + A^(rfit+THMF). (8) 



The first term is the cost of pre-calculating the HMF, 
with one iteration of constructing model visibilities and 
gridding to the uv plane times a computation overhead 
S associated with building a sparse matrix. The second 
term is the time required to perform the deconvolution, 
consisting of the time to apply the HMF thmf and the 
time to identify the next model component rgt times the 
number of model components N. The speed advantage 
of FHD comes from moving the costs involved in model- 
ing and gridding visibilities outside of the deconvolution 
loop. 

Table [1] lists the computational cost of each term in 
Equations [3 and [5] for our example MWA 32 data (fg)). 
The total time for deconvolving 20,000 components with 
the A-projection algorithm (ta) is more than two orders 
of magnitude longer than with Fast Holographic Decon- 
volution (tfhd)- 

These costs do depend on the details of the instrument. 
In general Tmodci and Tgrid will scale linearly with the 
number of visibilities and the size of the gridding kernel, 
while TfjMF scales linearly with uv coverage and with the 
square of the size of the gridding kernel. In practice there 
are many optimizations that can be applied to both al- 
gorithms, many of which will help both approaches. For 
our IDL implementation of Fast Holographic Deconvolu- 
tion we have been able to reduce the deconvolution time 
(2"'^ term of Equation [8|) to 23 minutes per polarization. 
Written in a more efficient programming language and 
parallelized, real-time deconvolution of the example in SjH 
should be achievable on a high-end desktop computer. 

4. RESULTS 

In this section we analyze a 5 minute snaps hot obser- 
vatio n by the 32 antenna MWA prototype (|Ord et al.l 
120101 ) to demonstrate the effectiveness of Fast Holo- 
graphic Deconvolution. Detailed descripti ons of the full 
as-bu ilt M WA facility can be f ound in (jTingav et al.l 
[201J) and ([Lonsdale et al.ll2009[ ). This observation is 
from 139 MHz-169 MHz with Pictor A at the center of 
the primary beam field of view (approx. 41° x 41° out to 
the first null) , and was taken on September 25*^, 2010 at 

Table 1 

Computational costs 



Single iteration 
Tmodci 123 CPU-s 
Tgrid 638 CPU-s 
Tflt 1.66 CPU-s 
THMF 3.21 CPU-s 



Full deconvolution loop 
iV 20,000 iterations 
S 9.59 
TA 4232 CPU-hours 

TFHD 29.1 CPU-hours 



5 



5:19 am local time (21:19 UT) on our 14*^ field expedi- 
tion. The Orion and Crab Nebulas are 40.5° and 67.9° to 
the north of field center in the first and second antenna 
side lobes, and the Galactic plane extends through the 
southern side lobes. The results of our FHD algorithm 
are shown in Figures [1] and [2] No attempt to deconvolve 
sources in the sidelobes of the array beam was made, so 
artifacts from the bright Crab nebula are visible in the 
Northern corners of panels (b)-(d). 

Panel (a) of Figure [1] shows the holographic dirty im- 
age prior to deconvolution, which is dominated by Pictor 
A and the array beam (PSF). We perform 20k iterations 
of deconvolution with a clean gain of 15% using only 
positive, centroided source components and removing all 
peaks above 600 mJy in the central region of the residual 
image. This peak detection threshold rises as sensitiv- 
ity decreases towards the edge of the field to maintain a 
constant minimum signal to noise. We then collect all 
components within 6 arc minutes to identify 2342 source 
candidates, shown in panel (b) [because each component 
is independently centroided, individual components of a 
real point source will differ slightly in location by small 
fractions of a pixel] . The residual image is shown in panel 
(c) on an expanded linear scale to highlight the fluctua- 
tions (no smoothing or filtering has been applied to any 
of these images) . There are very few features in the resid- 
ual image, the residuals are highly Gaussian with noise 
of a— 230 mJy in the center, and there is no evidence 
of the Pictor A array beam. As expected, the noise in- 
creases towards the beam edge where the sensitivity is 
down by a factor of 20. The flux density of Pictor A has 
been scaled to match the catalog value of 402 Jy to set 
an approximate overall scale. 

The source counts are shown in panel (a) of Figure [21 
The black solid line shows the counts from the 6'^ Cam- 
bridge survey at 151 MHz in the Northern hemisphere 
(|Hales et al.lll988D . overlaid by our source counts in red. 
Because we deconvolve sources out to near the first null 
of the tile beam, our sensitivity naturally falls towards 
the field edge. At the faint end of the distribution (be- 
low 3 Jy) we divide our source counts into rings of nearly 
equal sensitivity (dashed lines) . At the bright end (above 
3 Jy) we use the full field as the uncertainty becomes 
dominated by the low source counts. There appears to 
be good correspondence between the source counts from 
FHD and the expected counts, with the expected roll-offs 
in sensitivity in the outer edges of the field. We observe 
excellent agreement of all bright source candidates with 
the Molonglq Refe rence Catalog (MRC) of radio sources 
(|Large et al.lll98ll ). though we defer a thorough catalog 
comparison to future work with deeper integrations. 

HERA0 instruments such as the 32 antenna MWA pro- 
totype quickly reach the confusion limit (expected to be 
^300 mJy). This makes measuring the dynamic range 
of a deconvolution algorithm problematic using common 
measures such as the RMS in blank areas of the field. 
Within 7° of the field center, the source counts appear 
to follow the 6C catalog to 600 mJy This implies a 
minimum snapshot dynamic range for Fast Holographic 
Deconvolution of 800:1 in this data set. This dynamic 
range will naturally improve with synthesis rotation. The 
polarization performance is represented by the residual 

^ Hydrogen Epoch of Reionization Arrays [(reionization.org)] 



Stokes Q image in Figure [2Jd, with a fractional polar- 
ization of a few percent in the central region. With im- 
proved beam models and calibration we expect to achieve 
percent polarization residuals across the field. 

5. CONCLUSION 

A-projection and forward modeling deconvolution al- 
gorithms that use direction-dependent gain and polar- 
ization are very precise, but are much slower than the 
traditional CLEAN approach. In this paper we have in- 
troduced Fast Holographic Deconvolution (FHD), which 
uses the Holographic Mapping function to improve the 
speed of holographic deconvolution without sacrificing 
precision. We have demonstrated the power of the FHD 
algorithm using a 5 minute snapshot of data from the 32 
tile MWA prototype. 

In future work we plan to enhance the capabili- 
ties of our Fast Holographic Deconvolution algorithm; 
including parallelizing and optimizing our implemen- 
tation, deconvolving diffuse sources, and implement- 
ing w-projection to enable long widefield integrations 
(CornwcU et al.''2008'; using holographic formulation in 
[Morales fc Matcick 2009). Because w-projection de- 
creases the sparsencss of the Holographic Mapping func- 
tion (larger gridding kernels) it may be more efficient 
to do a joint deconvolution of a few intermediate length 
observations than one long track. We will also explore 
ways of using sparse matrix parallelization in addition to 
the more traditional frequency or time parallelizations. 
Finally, Fast Holographic Deconvolution lends itself nat- 
urally to modeling diffuse sources such as the galaxy, and 
we plan to explore multi-scale or other principle compo- 
nent deconvolution methods. 
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(a) Dirty Image (b) Sources 




(c) Residual Image (d) Restored Image 




Figure 1. Data from the MWA 32 antenna prototype analyzed with Fast Holographic Deconvolution. The data shown is purposefully a 
torture test for deconvolution algorithms — a 5 minute snapshot with Pictor A at the center of the ~ 41° X 41° primary beam FoV, the 
Orion and Crab Nebulas 40.5° and 67.9° to the north of field center (well off the image), and the Galactic plane extending through the 
antenna side lobes. Panel (a) shows the Stokes I dirty image, dominated by Pictor A and the snapshot PSF (array beam) shown on a log 
stretch to highlight low level detail. Deconvolution reveals 2342 source candidates in this 5 minute snapshot in panel (b) at upper right 
(log stretch). Panel (c) shows the residual image after deconvolution on an expanded linear scale to highlight the residual emission, and 
panel (d) shows the reconstructed image (panels b + c) on the same scale. There are false sources in the upper corners of the image duo to 
bright sources outside of the clean region which are not deconvolved. However, we see no residuals from Pictor A or other sources inside the 
deconvolved region. See the text for a detailed description of the features and Figure [2] for the log A'/ log 5 plot and the Stokes Q Image. 
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(a) Source counts (b) Stokes Q Residual Image 
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Figure 2. Panel (a) shows the logAf/logS plot for the sources in Figure[T] with error bars to indicate sample variance at high fluxes. We 
have divided the source statistics below 3 Jy into rings of increasing radius (dashed, dotted, dash-dotted). Careful investigation of each 
curve shows the turnover in the completeness increasing with radius as expected. In the central region of the image the sources follow the 
expected source counts until 600 mjy, which gives us a conservative snapshot dynamic range of >800:1. This dynamic range will naturally 
improve with synthesis rotation (dynamic range becomes hard to measure in longer MWA 32 observations due to the confusion limit). 
Panel (b) shows the residual Stokes Q image, with fractional polarization of a few percent in the central region (21° from antenna bore 
sight). With upgraded beam models and calibration we expect to achieve percent polarization across the field. 
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